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ABSTRACT 
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P^ ! The orientation of a massive binary undergoes a random walk due to grav- 

Q i itational encounters with field stars. The rotational diffusion coefficient for a 

circular-orbit binary is derived via scattering experiments. The binary is shown 



Q 



(N 
> 

\o 
in 
i> 

o 



o 



o 



to reorient itself by an angle of order ^JrajM during the time that its semi-major 
axis shrinks appreciably, where M is the binary mass and m the perturber mass. 
Implications for the orientations of rotating black holes are discussed. 



1. Introduction 



Two sorts of Brownian motion are defined by statistical physicists. Translational Brown- 
ian motion is the irregular motion exhibited by a massive particle as it collides with molecules 
in a fiuid. Collisions also affect the spin and orientation of particles; theories of rotational 
Brownian motion are concerned with the evolution of these quantities due to encounters, 
j^ \ Rotational Brownian motion was first discussed by Debye (1913, 1929) in the context of 

dielectric theory. The polarization of a dielectric material is a competition between torques 
^>< \ due to the imposed field, which tend to align polar molecules, and collisions, which tend to 

H \ destroy the alignment. 

The likely existence of black holes with masses much greater than those of stars has led 
to a renewed interest in Brownian motion among astrophysicists (Merritt 2001; Milosavljevic 
& Merritt 2001; Dorband, Hemsendorf & Merritt 2001). Translational Brownian motion of 
black holes in galactic nuclei may be directly observable (e.g. Reid et al. (1999); Backer 
& Sramek (1999)), allowing constraints to be placed on black hole masses. Translational 
Brownian motion may also affect the rate of tidal disruption of stars by a nuclear black 
hole (e.g. Young (1977)) or the efficiency with which a binary black hole can interact with 
stars as it wanders through a nucleus (Quinlan & Hernquist 1997; Milosavljevic & Merritt 
2001). In Paper I (Merritt 2001), the translational Brownian motion of a massive binary was 
discussed; it was shown that the amplitude of the wandering can be increased by a modest 
factor compared with that of a single point mass due to exchange of energy between field 
stars and the binary. 
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Here we consider the second kind of Brownian motion experienced by a binary, the 
irregular variation of the orientation of the binary's spin axis due to encounters. A single 
field star that passes within a distance ~ a of the binary, with a the binary's semi-major 
axis, will exchange orbital angular momentum with the binary, leading both to a change in 
the binary's orbital eccentricity as well as a change in the orientation of the binary's spin 
axis. Repeated encounters will cause the binary's orbital eccentricity to evolve (usually in 
the direction of increasing eccentricity) and will also cause the orientation of the binary to 
undergo a random walk. The first process has been discussed extensively (e.g. Mikkola & 
Valtonen (1992); Quinlan (1996)); the second, discussed here for the first time, will be called 
rotational Brownian motion, by analogy with the similar process that occurs in a fluid of 
polar molecules (e.g. McConnell (1980)). 

The orbital orientation of a black-hole binary influences the direction of the spin axis of 
the single black hole that forms after coalescence of the binary due to emission of gravitational 
radiation (e.g. Flanagan & Hughes (1998)). The spin axis, in turn, is thought to determine 
the plane of the inner accretion disk which forms around the black hole (Bardeen & Petterson 
1975), the direction of jets launched from the accretion disk (Rees 1978), etc. Thus changes 
in the orientation of a black-hole binary as it coalesces may have observable consequences. 

In §2 the conservation equations for interaction of a field star with a massive binary are 
presented; most of these have been seen before and were used in earlier studies to derive the 
rates of change of the binary's separation and eccentricity. The rotational diffusion equation 
is presented in §3 and the order-of- magnitude of the rotational diffusion coefficient (Ai9^) 
is derived. §4 describes the numerical derivation of {Ad"^) from scattering experiments. §5 
presents solutions of the diffusion equation and derives a relation between the degree of 
hardening of a binary and the change in its orientation. 



2. Encounter Kinematics 

Consider an encounter of a single particle (the "field star" ) with a binary. Assume that 
the binary remains bound during the encounter, appropriate if the mass of the field star is 
much less than the mass of each of the components of the binary or if the binary is sufficiently 
hard. At early times, the field star has velocity vq with respect to the center of mass of the 
field star-binary system and its impact parameter is p. A long time after the encounter, 
the velocity v of the field star attains a constant value. Conservation of linear momentum 
implies that the change (5V in the velocity of the binary's center of mass is given by 
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where ruf is the mass of the field star and M12 = Mi + M2 is the mass of the binary. This 
relation was used in Paper I in conjunction with scattering experiments to investigate the 
translational Brownian motion of a binary due to encounters. 

The energy of the field star-binary system, expressed in terms of pre-encounter quanti- 
ties, is 

1 2 1,. ..2 GM1M2 , , 

Eo = s'^/^o + 2^12^,' - ^^ (2a) 

1 A , ^/ \ 2 GM1M2 



with Vq = — (mj/Mi2)vo the initial velocity of the binary's center of mass and Oq the 
binary's initial semi-major axis. After the encounter, 

^1 / mf\ ^ GM1M2 



and E = Eq, so that 



a - GM1M2 V M12' ^ ' 



GM1M2 ^ ' 

where the latter relation assumes M12 ^ mj. This result was used by Hills (1983, 1992), 
Mikkola & Valtonen (1992), and Quinlan (1996) to compute the hardening rate, {d/dt){l/a), 
of a massive binary due to encounters. 

The angular momentum of the field star-binary system about its center of mass, ex- 
pressed in terms of pre-encounter quantities, is 

Co = mJl + ^jio + fx^2ho (5) 

where fj,i2 = M1M2/M12, io = P^o and ibo = 'C,bo/fii2 with Ci, the binary's spin angular 
momentum. Conservation of angular momentum during the encounter gives 

A4 = -^(l + ^)si (6a) 

/il2 V M12J 

^ -'^6£. (6b) 

Changes in \ii,\ correspond to changes in the binary's orbital eccentricity e via the relation 
e^ = 1 — £l/GMi2jj,'^a (Mikkola & Valtonen 1992; Quinlan 1996). Changes in the direction 
of ib correspond to changes in the orientation of the binary, leading to rotational diffusion. 
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3. Rotational Diffusion Equation 

Let F{6,(f),t)dQ be the probability that the spin axis of the binary is oriented within 
the solid angle dQ at time t. We seek an equation describing how F evolves with time 
due to encounters of field stars with the binary, if each encounter is assumed to produce a 
negligibly small change in 6 and 0. If we imagine that Cf, is initially oriented parallel to the 
^ = axis, and that the encounters leading to changes in the orientation of the binary are 
isotropic in velocity and direction, then F will evolve in such a way as to remain a function 
of 6 alone, F = F{6, t). These are the same assumptions made by Debye (1929) in his theory 
of the rotational Brownian motion of spherical molecules. Debye showed that the evolution 
equation for F, in the absence of an external torque, is 

dF I d f . AM^)dF\ 

^sm^V^— (7) 



dt sin edB\ A 89 

with (A-i?^) the rotational diffusion coefficient, 

{A^^)= f^{dn,dn')'d^dn'; (8) 

^(dfi, dQ')dQ' is the probability that, during a unit interval of time, a binary whose angular 
momentum Cb is directed toward dQ will reorient itself such that Cb lies within dQ', and i} 
is the angular separation between dQ and dQ'. 

From its definition, (A-i?^) is the sum, over a unit interval of time, of {6'd)'^ due to 
encounters with field stars. It can be expressed in terms of mean changes in the components 
of the binary's angular momentum as follows. Let the binary be initially oriented with its 
spin vector along the z axis, ibo = ibo^z- After an encounter, the angle i) between initial and 
final spins is 

cos^ = ^= ^^. (9) 

y H,x ^ '-b,y ^ '^b,z 

Expanding equation (9) to second order in i) and taking averages, we find 

(A.=) ^ (^fM±Ml ,10) 

which can further be expressed in terms of changes in the angular momentum components 
of the field star via equation (6b) as 

,^ „ m) {Ail) + {Ail) , . 






G'MiM2/ii2a(l - e2 



{{Ail) + (A^)] . (lib) 
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The last relations assume Mu ^ rrif. Specializing to the case of an equal- mass, circular- 
orbit binary, equation (lib) becomes 

(M^) = 16 —^ - — - — ^ — ^. 12 

Equation (12) is a prescription for computing the rotational diffusion coefficient from scat- 
tering experiments. This prescription will be applied in the next section. 

We can estimate the order of magnitude of the rotational diffusion coefficient by noting 
that ^JGM\2|a = Vun-, the relative velocity of Mi and M2, and that a single close encounter 
produces a change in l of the field star of order Vun x «• The rate of close encounters is 
~ 2nGMi2nfa/af with Uf and o"/ the number density and Id velocity dispersion of the field 
stars (Paper I). Thus (Af) ^ {aVun}'^ x 2nGMi2nfa/af and 

(A,.>«32.^^ (13) 

M12 (Jf 

where p/ = rnfUf. 

The same encounters that induce rotational diffusion will also cause the binary to 
harden, according to equation (4b). The hardening rate may be written 

d fl\ _ GpfH 



dt \aj af 

with H ^ 15 for a hard binary (e.g. Quinlan (1996)). Thus 



, . ^2v 327r rrif 
^ ' H M 



12 



d (\ 



a—- - 
at \a 



327r m 



f 



H M12 thard 



(15) 



with t^ard — ci{d/dt){l/a). Equation (15) suggests that rotational diffusion occurs on a time 
scale that is of order Mi2/mf times the hardening time, or that the binary will rotate by 
an angle of order ^JmjjM\2 in the time it takes to harden significantly. This result will be 
assessed more quantitatively in §5. 
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4. Numerical Computation of the Diffusion Coefficient 

The same set of scattering experiments described in Paper I were used here to evaluate 
{Ail) ^-nd (Aiy). As described more fully in that paper, field stars were treated as massless 
particles moving in the potential of the binary whose center of mass and orbital parameters 
remained fixed. Orbits were integrated using the routine D0P853 of Hairer (N0rsett & Wan- 
ner 1991), an 8(6)th order embedded Runge-Kutta integrator; each field star was assumed 
to begin at a position {x,y,z) = (oo,p, 0) and was advanced from r = oo to r = 50a along 
a Keplerian orbit about a point mass M12. The integrations were terminated when the star 
had moved a distance from the binary that was at least 100 times its initial distance with 
positive energy. The binary's orbit was assumed to be circular; the orientation of the bi- 
nary's orbital plane with respect to the initial velocity vector of the field star, and the initial 
phase of the binary, were chosen randomly for each integration. Unless otherwise indicated, 
distances and velocities are given below in program units of a and {GMi2/ay^'^ respectively. 
The orbit integrations, and the expressions given above relating changes in the field-star's 
angular momentum to the orientation of the binary, were checked by carrying out a lim- 
ited set of experiments using a fully general, three-body integrator and various values for 
mf/Mi2. The expressions derived above for Si) etc. were found to be accurately reproduced 
in the limit of small mf/Mi2. 

Consider first the rotational diffusion coefficient corresponding to interactions with stars 
of a single velocity V at infinity. Multiplying {6'dy in a single encounter by 27ipnfV, with p 
the impact parameter and n/ the density of field stars, and integrating over p gives 

{6^\V,p^a.)) = 2nnfV dp pW (16a) 



rPmax I ^p _|_ (5£2 \ 

= 2TmfV / dpp { ''" ''^ (16b) 

(cf equation 12), where 511,51'^ are defined as mean square changes at fixed (V,Pmax), 
averaged over many encounters in which the orientation of V with respect to the binary's 
spin axis is taken to be random. Figure 1 shows the distribution of field-star velocity changes 
as a function of P ioi V = O.SVbm- 

We define a dimensionless, velocity-dependent diffusion coefficient D^ as 

{5'd'^(y,Pmax)) = ™^«/Vbin f TJ^ J ^ij ( V, Pmaa; ) , (17a) 
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Fig. 1. — Distribution of field star angular momentum changes for scattering experiments 
with V/Vbin = 0.5. Units are G = M12 = a = 1. Each plot corresponds to 5 x 10"^ scattering 
experiments within some range of impact parameters [pi,P2] in units of a. (a) [4,6]; (b) 
[2,4]; (c) [1,2]; (d) [0.6,1]. 



y \ rPmax/a 2 [Sil + Si"^ 



DR{V,Praa.) = 16 -— / d ('-) ^ ' (17b) 



with Vun = \jGMi2/a. Figure 2 plots D^ as a function of V and Pmax- The most ef- 
fective encounters are those with small V and Pmax'-, encounters with impact parameters 
p > ^/2GMi2a/V have a closest-approach distance to the binary oi r-p > a and produce 
almost no change in the binary's orientation. Hence Dr reaches a maximum value at 
p K, v^2GMi2a/V" before levelling off at larger p. 

Multiplying {dd"^) by ff{vf), the distribution of field-star velocities, and integrating dvj 
gives (Ai9^). We assume a Maxwellian distribution of field star velocities, 



ffi^f) 



(27ra2)3/2 
so that 



-vj/2aj 



{M') = -^J rfV\/V^^/-?(5^^). (19) 

(ivro-^) ' Jo 



This can be written in terms of a dimensionless coefficient L, 



2, _ ruf Gpftt 



(A^^) = j^^^^^L, (20a) 

Mi2 (Tf 

fOO 

(20b) 



L = V2^S'^ dz z'^e-''/^DR{zVun/S,Rpf). 
Jo 



L is defined in analogy with the dimensionless coefficients {H, J, K) that describe respec- 
tively the rates of change of the binary's energy, the rate of mass ejection by the binary, and 
the binary's eccentricity growth rate (cf. Quinlan (1996)). Note that L is a function of two 
parameters, L = L{R, S), where 

J3, _ Pmax PmaxO'f _ Vhin . ^ 

Pf GMi2 CTf 

R is the maximum impact parameter in units of the radius of gravitational influence of the 
binary and S is the dimensionless hardness of the binary. Figure 3 plots L{R, S). At a given 
hardness S, L reaches its limiting value Loo = L{oo, S) hj R ^ 3/S. For all but the softest 
binaries, this limiting value is reached before R = 1, which is roughly the value of R we 
expect for a massive binary at the center of a steeply- falling density cusp (Paper I). Hence 
we ignore the dependence of L on Pmax in what follows and write L{S) = L{oo,S). L{S) 
varies from ~ 25 for the softest binaries {S ^ 1) to ~ 60 for S ~ 10. 
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Fig. 2. — Velocity-dependent rotational diffusion coefficient Dji(y,Pmax) defined in equa- 
tion (17b). Pmax is the maximum impact parameter for encounters and a is the bi- 
nary semi-major axis. Different curves correspond to relative velocities at infinity of 
V/Vbin = 0.05,0.1,0.2,0.3,0.5,1,2,10; the lowest (highest) V produces the highest (low- 
est) Dji at large Pmax- Points are averages computed from the numerical integrations; curves 
are spline fits. 
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Fig. 3. — Rotational diffusion coefficient L (equation 20) for a circular-orbit binary. Curves 
are marked by Vbin/c/ = 1, 2, 4, 6, 10. 
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5. Solutions of the Evolution Equation 

Equation (7) describes the evolution of F{9,t), the probability distribution for the bi- 
nary's orientation, given the rotational diffusion coefficient (Ai9^). The latter (equation 20) 
is a function of the binary separation a, directly via the factor Gpfa/af, and indirectly via 
the weak dependence of L on binary hardness (Figure 3). Furthermore the timescale for the 
binary to reorient itself is shorter than the hardening time \{l/a){da/dt)\^ (cf. equation 15), 
hence the change in a can not be neglected when solving for the evolution of the orientation. 

We nevertheless begin by considering an idealized model in which (A-i?^), a, and the pa- 
rameters {pf, o"/) that describe the background are assumed constant in time. The evolution 
equation (7) becomes 



dF 

97 



28^ 



(1-/^^) 



dF 
dp 



where p = cos 6, t = t/to, and 



M 



12 



^/ 



(22) 



(23) 



ruf J GpfaL 
Consider the evolution of 'p{t) = 27r J-^ F{p' , t)p'dp', the expectation value of p. We have 



dp 

dr 



"1 Qp 

2n I -;:—p'dp' 
'_! or 



n 



' ,d 

p 



1 



dp' _ 



^^-^'^m 



-27r 



p'Fdp' 



-p 



(24a) 
(24b) 
(24c) 



or 



p{t) = PqC \ (25) 

Hence to is the time constant for relaxation to a uniform distribution of orientations, /^ = 0, 
given a fixed binary separation and unchanging stellar background. It may be expressed in 
physical units as 



to = 6.1 X 10 V 



M 



12 



106 



mf 



L 

50 



-1 



Pf 



lO^MeMpc" 



-1 



^/ 



200 km s^ 



1 pc 



-1 



(26) 



Next consider the more interesting case where a, as well as the parameters pf and o"/ 
that describe the field stars, may be changing with time. This time dependence is not easily 
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specified but we can make progress by changing evolution variables from t to a; = log(a/ao). 
First combining equations (7) and (20a), we find 



Now changing variables, 



(cf. equation (14)) and so 



dF 

dt 


L nif Gpfa d 
4 Mi2 af d/x 


'(1 


2^dF 


- 


dF 

dt 


dF dx dF 1 da 
dx dt dx a dt 


dF GpfHa 

dx af 


dF 
dx 


1 L ruf d 
AHMi2dii _ 


(1- 







(27) 



(2^ 



(29) 



This equation gives the evolution of the binary's orientation in terms of changes in its 
semi-major axis a, with no explicit dependence on the parameters that describe the stellar 
background. Ignoring the weak dependence of H and L on a, appropriate for a hard binary, 
we find as before 

'' a 



p{a) = Pq 



ao 



1 L "»/ 



(30) 



The exponent in this expression is of order nif /M12 ^ 1. Hence we can write 

1 L ITif 



log 
2HMi2 \ao 



(31) 



where Pq has been set to unity, corresponding to an initial orientation parallel to the ^ = 
axis. 

For the case of binary supermassive black holes, if we define Oq as the separation when 
the binary first forms a bound pair, we expect gravitational radiation coalescence to occur 
when a/oo ~ 10~^ (Merritt 2000). Since L/2H ^ 2 for a hard binary, we find for the 
expectation value of p at coalescence: 



71^ 1 -10 



riy_ 



(32) 



12 



Writing 59 
becomes 



2(1 — p), the rms change in the angle defined by the binary's spin axis, this 
5^~W=^^, m/<Mi2. (33) 



Ml 



12 
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For iTif /Mi2 = 10 ^, e.g. a binary of mass lO^M© surrounded by IMq stars, this predicts 
66 ~ 0.25°. For a binary of so-called "intermediate-mass" black holes with M12 ~ IO^Mq, 
59 ^ 8°. 



Milosavljevic & Merritt (2001) give Information about the rotational Brownian motion 
of a massive binary in a set of A^-body simulations. They show (in their Figure 11) the 
evolution of the orientation of an equal-mass, circular-orbit binary in three simulations with 
mf/Mi2 = {164,328,655}. The simulations all follow the evolution of the binary until 
a/ao ~ 10~^; hence equation (32) predicts 69 ^ {14°, 10°, 7°}. These values are in excellent 
agreement with the angular deflections shown in Figure 11 of Milosavljevic & Merritt (2001). 

Under what circumstances could the roeorientation angle of a supermassive black hole 

1/2 
binary be much greater than the modest values predicted here? Since 69 oc mJ , one 

possibility is for the perturbing objects to be much more massive than IMq; for instance, they 

could be supergiant stars (m/ ~ lO^M©), star clusters {rrif ^ IO^Mq), or giant molecular 

clouds {rrif ~ 10'^Mq). Such physically large objects would presumably produce most of 

their perturbing torques when the binary separation was greater than or comparable to their 

own size. 

An even more extreme possibility is that the perturbing objects are themselves super- 
massive black holes. Consider for example a binary supermassive black hole at the center 
of a giant elliptical galaxy. The stellar densities at the centers of such galaxies are low and 
one might expect the coalescence of the binary to stall due to a shortage of interacting stars 
(e.g. Milosavljevic & Merritt (2001)). In this case, further decay of the binary might await 
the infall of a third and subsequent black holes from later galactic mergers (e.g. Valtonen 
(1996)). Mass ratios irtf /M12 would then be of order 10~^ — 10° and signficant reorientation 
of the binary could occur before coalescence. 

Orientations of radio jets in Seyfert galaxies are known to be almost random with respect 
to the plane of the stellar disk (e.g. Kinney et al. (2000)). If jets are launched parallel to 
the spin axes of rotating black holes, this implies that the spins of nuclear black holes 
bear almost no relation to disk spins. While the modest reorientation rates derived here are 
probably insufficient to explain this phenomenon, a related mechanism may suffice. The spin 
orientation of a supermassive black hole that forms via repeated mergers with smaller black 
holes would undergo a random walk due to the random orbital inclinations of the infalling 
black holes. If nuclear black holes in Seyfert galaxies grew through multiple mergers before 
the formation of the stellar disks, one might expect their spin axes to be uncorrelated with 
those of the stars and gas. 
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